Electromagnetically induced transparency in a superconducting three-level system 
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When a three-level quantum system is irradiated by an intense coupling field resonant with two 
of the three possible transitions, the resonant absorption of the system from its ground state by 
an additional radiation field is suppressed. This effect, where the population is trapped in the 
ground state, is known in quantum optics as "electromagnetically induced transparency". When 
the coupling field is detuned from resonance, the resonant absorption peak splits to form an " Autler- 
Townes doublet". We observe these phenomena in a superconducting Josephson phase qubit, which 
can be considered an "artificial atom" with a multilevel quantum structure. These observations are 
qualitatively described by a simple model restricted to three energy levels. A full solution of the 
master equation including higher levels provides excellent agreement with all the experimental data. 
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Superconducting qubits [l| have recently been em- 
ployed as a testing ground for quantum mechanics in a 
nearly macroscopic system. The interaction of the ef- 
fective two-level system with a resonant cavity has at- 
tracted a lot of attention 0, 0, 0, @] • However, Josephson 
junction-based quantum systems can also possess higher 
quantum states or energy levels, beyond the two- level ap- 
proximation, and these have received much less attention. 
For driven Rabi oscillations between the lowest two levels 
of a phase qubit, these higher states can reduce the ap- 
parent Rabi frequency at large drive powers @, 0] as the 
energy leaks out of the subspace spanned by the ground 
state |0) and the first excited state |1). Recently, the sec- 
ond excited state |2) has been used in order to help char- 
acterize the fidelity of single quantum bit (qubit) gate 
operations 0. 

In atomic physics and quantum optics, various exper- 
iments have utilized the naturally occurring multilevel 
state structure of "real" atoms. A system accessing just 
three energy levels can display a rich variety of phenom- 
ena. Coherent population trapping, electromagnetically 
induced transparency (EIT) [lCj , Autler-Townes splitting 
TO], and stimulated Raman adiabatic passage [12|], have 



been thoroughly investigated with atoms, and recently, 
also with solid-state quantum dots is in- 

triguing to demonstrate some of these phenomena using 
macroscopic quantum states in a superconducting-based 
framework. Phase qubits in particular are well-suited 
for this purpose because they have a simple ladder-type 
energy level structure as well as a measurement process 
that is fast and state specific. 

In a phase qubit (see Fig. 1(a)), a single Josephson 
junction, capacitively shunted with Cj, having an effec- 
tive linear inductance Lj = ($o/27r) /Ej, where Ej 
is the Josephson energy, and <&o is the flux quantum, 
has been inserted into a superconducting loop. The 
loop inductance L > Lj is chosen so that local min- 



ima are formed in a one-dimensional energy potential 
E pot = ($ - $ 0Xt ) 2 /2 J L - £:jcos(27T$/$ ) controllable 
by an externally applied flux $ ex t. For this experiment, 
there are approximately ten energy levels residing within 
one of these local minima as shown in Fig. 1 (b) . 

Consider the three lowest energy levels in Fig. 1(b). We 
denote the transition frequencies between energy levels i 
and j as Uji — 2wfji. Two microwave fields are present; 
a weak probe tone lu p that is nearly resonant with u>iq 
and detuned an amount A p = luiq — u) p , and a strong 
coupling tone lu c that is nearly resonant with u>2i (de- 
tuned by A c = W21 — w c ). The amplitude of these tones 
leads to corresponding Rabi "flopping" frequencies de- 
noted by and Q\q. In atomic physics, the occupation 
of higher energy levels will often decay into a continuum 
of states, but here, the occupation of the higher levels 
will decay merely to the lower levels of this ladder or 
cascade- type configuration as seen in Fig. 1(b). This can 
be contrasted with the so-called "lambda" configuration 
familiar in atomic systems, where the state |2) would be 
lower in energy than the state |1) . We denote the inter- 
level relaxation rates in our system by T2i and Tio- The 
rate is a forbidden dipole transition due to the an- 
harmonicity of the energy potential. 

In the qubit 's eigenbasis |n), the Hamiltonian takes 
the form H = ^ n=Q E n \n)(n\ + fl p cos(u> p t)\l) (0\ + 
f2 c cos(a; c t)|2)(l|+ h.c. In the rotating frame, neglecting 
counter-rotating terms, the Hamiltonian takes the form 



H 



o n p /2 
o n c /2 



o 



(1) 



This Hamiltonian is familiar in atomic physics and has 



been used to describe coherent population trapping [15 1 
and EIT [l(| • Analogous phenomena have been predicted 
for superconducting qubits [3,0], however, experimen- 
tal measurements have not been forthcoming until now 
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FIG. 1: (a) Schematics of thephase qubit and a nearby dc- 
squid for single-shot readout [3, [ly] and on-chip flux coil; (b) 
Energy levels in the local potential minimum at $ ~ $0 given 
for the flux bias point $ C xt ~ —0.42 $0 used in the experi- 
ment. The zoom illustrates the cascade configuration of the 
three lowest levels used to generate electromagnetically in- 
duced transparency and the Autler-Townes splitting. The 
(strong) coupling field u c connects connects energy levels |2) 
and |1), and the weaker probe field f p , whose resonant ab- 
sorption becomes blocked due to the coupling field, connects 
|1) and |0). The corresponding Rabi drive amplitudes are de- 
noted by Q c and fi p , respectively. Generally, both tones are 
detuned from the respective transitions by A c and A p ; (c) 
Two-tone microwave spectroscopy to reveal the two lowest 
transitions /10 = 8.135 GHz and /21 = 7.975 GHz (see text 
for details). 



.18]. 

When both the coupling and probe tones are resonant 
with their corresponding transitions (A c = A p = 0), 
one of the eigenstates of Eq. |T]) is the "dark state" 
*o = cos(6)|0) - sin(6)|2), where tan(9) = n p /n c . 
Note that this expression does not include the interme- 
diate state |1). Thus, under these conditions, the system 
can no longer absorb radiation at the probe tone and 
becomes "transparent". Next, we describe experiments 
demonstrating this phenomenon for the multilevel super- 
conducting phase qubit shown in Fig. 1. 

The qubit sample was made by the use of optical 
lithography and standard A1-A10 X -A1 ion-mill tunnel 
junctions on a sapphire substrate. As compared to 
Ref. 0], the junction size was reduced by a factor of two 
in order to reduce the number of microscopic TLS defects 
within the tunnel barrier [l9| . The experiments were per- 



formed at 25 mK in a dilution cryostat. The qubit state is 
measured with a nanosecond- wide flux pulse as described 
in Ref. [H . The measure pulse reduces the potential bar- 
rier AU (Fig. 1(b)) so that the resultant energy level at 
the top of this barrier, if occupied, has a high probability 
of tunneling out of the local minimum. Calibration of the 
measure pulse amplitude for the first excited state gives 
its population Pi , and also leads to a residual tunneling 
probability Pq out of the ground state of approximately 
5%. Using this calibration, we also measure the popu- 
lation Pi + P2, when the second excited state |2) is oc- 
cupied. To measure Pi independently, we calibrate the 
measure pulse amplitude so that the residual tunneling 
probability Pi out of |1) is approximately 5%. 
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FIG. 2: Driven Rabi oscillations between energy levels |2) 
and |1). (a) Excited state probability P2 as a function of the 
width of a pulse resonant with /21 (this pulse follows a n pulse 
resonant with fio, see text); (b) the frequency of the Rabi 
oscillations Q21 as a function of the coupling drive amplitude 
Q c . The dashed line is a fit to the linear portion of the data. 



We operate the phase qubit at a center frequency 
fio = 8.135 GHz where there are no visible TLS defects 
over a 1 GHz bandwidth. At high microwave powers, we 
spectroscopically observe a two-photon absorption peak, 
/20/2 = 8.06 GHz. In order to identify the transition fre- 
quency /21 , we simultaneously apply two weak microwave 
tones. We set the probe tone f p to resonantly pump the 
/xo transition while we vary the coupling tone f c . We 
obtain a peak in the spectrum when f c = /21 = 7.975 
GHz for different measure pulse calibrations as shown in 
Fig. 1(c). The full spectroscopic data (not shown) for 
fio and /21 as a function of the dc flux bias can be fit 
well to theory in order to extract the qubit parameters 
Lj and Cj, including any flux offsets. We also perform 
time-domain measurements (Fig. 2) by first applying a 
7r pulse at fio in order to populate |1) followed immedi- 
ately by a pulse resonant with /21 of varying length. We 
obtain coherent Rabi oscillations of the second excited 
state P2, where the population flops between |2) and |1). 
In Fig. 2(b), we show the Rabi frequency O21 of these os- 
cillations for several different coupling drive amplitudes 
£l c . A linear relationship between the Rabi frequency and 
drive amplitude survives up to about ^21 ~ 70 MHz, in- 
dicating that a three-level model is sufficient to describe 
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the dynamics up to this drive 0, 0] ■ 

We made independent measurements of the inter-level 
decay rates F 2 i and Tio (Fig. 1(b)) by preparing popu- 
lation in either |2) or |1), using sufficiently slow (~ 10 
ns wide) 7r pulses in order to avoid population leakage 
to the higher levels, and observing energy decay in time 
domain. This gave T 2 i = 11 MHz and r i0 = 7 MHz. 
The widths <5/;j of the spectroscopy peaks in Fig. Q] (c) 
give the pure dephasing rates = 3.5 MHz, r£, = 6 
MHz, according to 
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FIG. 3: Demonstration of electromagnetically induced trans- 
parency (EIT) in a multilevel phase qubit system. Plotted is 
the excited-level population Pi + P2 at different amplitudes 
of the coupling field, which cancels the resonant absorption 
of a weak probe field Q p = 3 MHz. (a) experiment; (b) the- 
ory. Solid lines are the full simulation (see text), dashed lines 
(shifted vertically) represent the analytical line-shape from 
Eq. (2); (c) splitting as a function of the independently mea- 
sured Rabi frequency f22i (from Fig. 2(b)). The solid line 
demonstrates the expected identity between the two quanti- 
ties. 

As discussed earlier, the "dark state" is expected to oc- 
cur when both the coupling and probe tones are resonant 
with their corresponding transitions (A c = A p = 0). We 
expect this condition to be visible as a reduced popu- 
lation of the first excited state Pi. In our cascade level 
configuration, energy decay from |2) reduces the effective- 
ness of the transparency by populating |1). Therefore, we 
increase the coupling tone amplitude so that f2 c 3> Q p , 
and the dark state approximates the ground state. We 
find that the probe tone absorption peak begins to split, 



forming a dip at f p — fiQ. The fact that the linewidth 
of this dip is broader than the linewidth of the /10 ab- 
sorption peak, is another indication of the middle level 
|1) receiving population. 

The characteristic features of EIT can be investi- 
gated further by determining the stationary states of the 
driven multilevel system by use of the Markovian mas- 
ter equation for the qubit's reduced density matrix p: 
ihp = [H,p] + C[p\. The relaxation term C[p] includes 
the measured inter-level relaxation rates T21 and Tio as 
well as pure dephasing. If we restrict ourselves to the 
first three energy levels, neglect fast oscillating terms, 
assume A c = 0, a relatively large coupling amplitude 
(fl c 3> (ftp, T 2 i, Tio)), and solve the resulting algebraic 
equations, we obtain a measured excited state popula- 
tion: 



Pi 



n 2 

P-2 = Pll + P22 = X 
JL 10 



4A£(A 



10 



rr ) + ^(Fi 



2A 2 



(2) 



16A4- 



[f^ + 4A 2 iAi 



8AM- 



-2 (All 



Afo)] 



where Ay = ry+I*. The splitting of the absorption peak 
has a characteristic line shape which is not just a sum of 
two Lorenzians. The stationary off-diagonal elements of 



the density matrix when A c 



= indicate the co- 
poi = 2iTiopii/n p . 



herence of the trapped population 
In Fig. 3(b), we find good agreement between the exper- 
imental data and the predictions from Eq. (2). In order 
to get the best fit, we have adjusted the value of Q c by 
about +10% from that of the independently measured 
Rabi frequency f2 2 i. 

The size of the splitting can be understood in terms of 
the Autler-Townes doublet 11, 21, 221 known in atomic 



physics. If we treat the amplitude of the probe tone as 
a small perturbation in Eq. ([1]), we obtain eigenenergies 
of the driven system: e± = wio + A c /2 ± y / A^+Of/2. 
These energy levels are excited from the ground state 
using the probe tone. When A c = 0, the doublet is 
spaced by the coupling Rabi amplitude fl c . In Fig. 3(c), 
this splitting is plotted as a function of the independently 
measured J7 2 i (from Fig. 2 (b)), confirming their identity. 
By scanning the frequencies of the probe and coupling 
tones, one expects to find that the Autler-Townes doublet 
will exhibit an avoided crossing centered at (u>xo, w 2 i). 
We display the results from this type of a measurement 
in Fig. [4] with the Autler-Townes eigenvalues plotted as 
white solid lines. Here, the experiment is performed with 
the same coupling amplitude as for the middle (black) 
curve shown in Fig. 2(a). Again we find good agreement 
for the simple three-level model when operating in the 
linear regime of the driven Rabi oscillations at £l 2 i = 36 
MHz (Fig. 2(b)). 

The vanishing intensity of the lower right leg of the 
Autler-Townes avoided crossing in Fig. [5] can be ex- 
plained by considering the matrix element sin((/>/2) con- 
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FIG. 4: Autler-Townes splitting in the three-level phase qubit, 
for a corresponding coupling field amplitude of tt c — ► Q.21 ~ 
36 MHz, and probe amplitude Q p — 5 MHz. Resonances in 
the driven system exhibit an avoided crossing as a function 
of the detuning of the coupling field (vertical axis) and the 
probe field (horizontal axis). The white lines mark the ex- 
pected transitions from the three-level Autler-Townes model, 
(a) measured excited-level population Pi + P2', (b) five-level 
simulation. 



necting the driven eigenstates with the ground state, 
where the mixing angle is given by tan(</>) = Q C /2A C . 
As the detuning A c increases, this matrix element de- 
creases so that the intensity of the split peak on the 
right also decreases, in agreement with the experimen- 
tal data. However, the opposite, upper left leg shows a 
decreasing then increasing intensity with detuning. This 
can be attributed to an additional two-photon process, 
which becomes favorable when the probe tone reaches 
the two-photon resonance at 8.06 GHz, also causing an 
enhancement of the base level. 

For the Hamiltonian of Eq. (JTJ, we have assumed that 
both microwave fields couple only to their intended tran- 
sitions, namely the coupling tone uj c drives W21 and the 
probe ujp drives u>io ■ In reality, the microwave fields cross- 
couple to both transitions. The strong coupling tone 
drives the f±o transition with a cross-coupling strength 
Q®, which contributes to the occupation of the first ex- 
cited state, Pi = |Oji/ (A c + ojio — cj 2 i) 2 + f^c > ris- 
ing the base-level in Fig. 3(b) when far from /10. For our 
weakly anharmonic potential, we estimate that the cross- 
coupling is Vl x c ~ Q c /V2. When Q 2 i = 36 MHz, (mid- 
dle (black) curve in Fig. 3, Fig. 4), the cross-coupling is 



negligible. However, for strong coupling amplitudes cor- 
responding to £l c = 66 MHz (top (red) curves in Fig. 3) 
approaching the limitations of the thrcc-lcvcl model, we 
must include the cross-coupling in order to correctly 
account for the position and symmetry of the Autler- 
Townes doublet. 

Thus far, our simple three-level model has qualitatively 
accounted for the EIT behavior and the resulting Autler- 
Townes splitting in good agreement with the experimen- 
tal data. Next, we provide a more accurate description 
of the system based on a full simulation that can quan- 
titatively account for the multilevel nature of the phase 
qubit system. We solve the Markovian master equation 
including the first five energy levels of the system and 
also including the cross-coupling of the drive tones. The 
results are shown along side the experimental data in 
Fig. 3(b) and Fig. 4(b). We find excellent agreement 
without any fitting parameters. The simulation correctly 
captures the asymmetry of the splitting (Fig. 3(a)), which 
is a combination of the effects of the higher levels and 
cross-coupling, as well as the intensity vanishing of the 
left and right branches as a function of detuning as seen 
in Fig. 4(b). 

In summary, we have observed phenomena character- 
istic of three-level systems familiar in atomic physics, 
namely electromagnetically induced transparency and 
Autler-Townes splitting, but in a superconducting 
Josephson junction-based quantum system. The results 
contribute to a general scientific effort that seeks to 
demonstrate quantum mechanical behavior in progres- 
sively more macroscopic and diverse systems. They also 
pave the way towards quantum information processing 
using higher-dimensional Hilbert spaces [23| . 
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